Operator Space Entanglement Entropy in XY Spin Chains 
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The complexity of representation of operators in quantum mechanics can be characterized by the 
operator space entanglement entropy (OSEE). We show that in the homogeneous Heisenberg XY 
spin 1/2 chains the OSEE for initial local operators grows at most logarithmically with time. The 
prefactor in front of the logarithm generally depends only on the number of stationary points of 
the quasi-particle dispersion relation and for the XY model changes from 1/3 to 2/3 exactly at the 
point of quantum phase transition to long-range magnetic correlations in the non-equilibrium steady 
state. In addition, we show that the presence of a small disorder triggers a saturation of the OSEE. 
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Complexity of many-body quantum systems has two 
remarkable manifestations: on one hand, certain strongly 
correlated (entangled) many-body states (e.g. cluster 
states [H) can be used to perform universal quantum 
computation and thus yield exponential gain over best 
known classical algorithms for certain tasks. On the 
other hand, weak entanglement (finite range of quantum 
correlation) is necessary for quantum states in order to 
be efficiently described classically 0. The latter fact is 
the reason behind the success of numerical methods, such 
as density matrix renormalization group (DMRG) 
which approximate quantum states by their most entan- 
gled components, e.g. as matrix product states (MPS) 
where the number of components depends on the quan- 
tum entanglement of the state. While the ground states 
of one-dimensional many-body systems are weakly en- 
tangled [1, d, 01 and can be efficiently computed by the 
DMRG, the time evolution of generic quantum states 
produces entanglement which results in a growing need 
of resources (see e.g. 0)- Such algorithms are only effi- 
cient if the computational costs grow at most polynomi- 
ally with time or, equivalently, the entanglement entropy 
which is a measure of quantum entanglement, grows no 
faster than logarithmically - which is the case only for 
very particular initial states and model systems. 

Another option is to simulate time-dependent oper- 
ators in Heisenberg picture using the same principles 
where the operators are represented as products of ma- 
trices, i.e. matrix product operators (MPO) |8|. This 
approach can be efficient even when the simulation of 
(pure) quantum states is not as is the case in transverse 
Ising spin chain 0, 0, OH 11 1 . In this context the relevant 
quantity determining the efficiency is the operator space 
entanglement entropy (OSEE) [12J which is an operator- 
space analogue of the usual bipartite entanglement en- 
tropy and presents a complementary quantification to 
other measures of entanglement in operator space 
As in the evolution of quantum states, similar limita- 
tions also apply to operators which accounts to the inef- 
ficiency of time-evolution of generic operators attributed 
by the exponential growth of resources in time. Never- 



theless, it was shown 0, [H| that local operators - on 
the contrary - can always be simulated efficiently in the 
integrable transverse Ising model with the OSEE being 
either finite (increasing logarithmically) for initial oper- 
ators which are represented repectively as products of 
a finite (infinite), number of Majorana fermions [having 
finite (infinite) index]. 

In this letter we shall consider OSEE of local Heisen- 
berg operators of infinite index in the quantum XY spin- 
1/2 chain, or any translationally invariant spin chain 
which can be solved by Wigner- Jordan transformation. 
We show that OSEE in such models generally increases 
logarithmically in time, where the prefactor is given uni- 
versally as 1/6 times the number of stationary points of 
the quasi-particle dispersion relation. We identify two 
generic regimes in the XY model with prefactors 1/3 and 
2/3 exactly corresponing to two quantum phases of a non- 
equilibrium steady state of open XY spin chain far from 
equilibrium [l5| and the regime of the gapless XY model 
with a prefactor 1/6. Investigating the effect of disorder 
we find, as an interesting consequence of Anderson-like 
localization in operator space, that OSEE saturates in 
time even for the infinite index initial operators. 

Dynamics of the quantum XY spin 1/2 chain of length 



n is described in terms of Pauli operators o ^ 
{1, 2, . . . , n} by the Hamiltonian 
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which is conveniently expressed as a quadratic form H 
w ■ Hw in terms of 2n Hcrmitian Majorana operators 



W 2 j-1 



(n 

l<3 



IP 

l<3 



(2) 



obeying the anticommutation relation {"Wj, w{\ = 26 ji, 
and 2n x 2n antisymmetric Hermitian matrix H. For 
the XY-model ((TJ the only upper-diagonal elements of 
H read H 2j>2j+1 = -(i/2)±±2 H 2j -i,2j +2 = (i/2) V» 
H 2 j-i, 2j = -\i/2)h for j e {1, . . . , n}. 
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As shown in 12j, we identify 4™ dimensional Pauli 
algebra with a Fock space of operators describing 2n ad- 
joint fermions (a-fermions) , with an orthonormal canon- 
ical basis \Pa) t = \w^w% 2 ■■■w%* n ) t , a 3 £ {0,1}. We 



define a set of adjoint annihilation linear maps Wj 



de- 



fined as Wj\P a )^ = aj\wjP a )^ which satisfy the canonical 

anticommutation relations {u^, — 0, {wp wj} = 5ji 
for j, I £ {1, . . . , 2n}. The Heisenberg dynamics in the 
adjoint (operator) space is given by a formal Schrodinger 
equation i(d/dt)A(t) = HA(t), or for the Majorana gen- 
erators i(d/dt)w(t) = 4Hw, which defines the adjoint 
hamiltonian H = —&dH = [»,H] = —Aw} ■ Hw. 

The operator space entanglement entropy (OSEE) of 
an arbitrary operator A = ^ Q a a P a is defined as the 
bipartite entanglement entropy of the adjoint state (a- 
state) \A)q — J2 a a a\Pa)^ in the operator space. The 
transformation between the physical basis {ctJ 1 • • •o£ n } 
and the basis {P a } is a simple permutation (with mul- 
tiplications by ±i or ±1) and despite the nonlocality it 
maps the first n/2 spin operators to the first n a-fermions 
and vice versa. Therefore, the OSEE of initially product 
operator A can be calculated by essentially following Ref. 
[H using the correlation matrix Tji = ^(AlwjW^A)^ for 
j, I £ {l,...,n}. The correlation matrix is Hermitian 
and the OSEE is calculated from its eigenvalues jj as 



S = J2H 2 ( 7j ) 



(3) 



3=1 



for H2(x) — —a; In a; — (1 — x) ln(l — x). 
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FIG. 1: Schematic examples of finite (i) and infinite (ii) index 
operators in operator Fock space. 

In this work we will only be interested in the time evo- 
lution of initially local operators i.e. products of a finite 




FIG. 2: OSEE of \F(t))* in the XY(7,/i) model with open 
boundaries at t — 800 and n = 1800. The dashed line denotes 
the critical line h c = |1 — j 2 \. 



the simplest infinite-index operator F corresponding to a 
half- filled Fermi sea of a-fermions (Fig. [1] case ii). Despite 
the nonlocality of the operator F — \ n / 2 a\ ■ ■ ■o~^, 2 , its 
entanglement properties are similar to those of a local 
operator a*, 2 = ]T n l 2+1 Fw n . 

The correlation matrix T for a time dependent ini- 
tially product a-state is calculated employing 
the Heisenberg picture in the operator space using 



t{A(t)\w} Wl \A(t)) t = t {A\w}(t) Wl (t)\A) t where w 3 {t) 



obtained from (d/dt)wj 
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H] as 



w{t) — $ T w where 



-4UH 
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Note that the matrix $ is real for any w-quadratic Hamil- 
tonian with Hermitian anti-symmetric matrix H. The 
correlation matrix used to calculate the OSEE ([3]) thus 
reads I>(t) = Y%Li %3^Pl*{ A H™p\A) r 

In Fig. [2 showing phase diagram 7 — h of OSEE in 
the XY-model at fixed large time t = 800, two distinct 
regimes are identified (plus a critical regime at h = 1). 
We will later show that these two regions are bordered 



number of Pauli operators a'J . This implies that, either b y the critical field strength h c = \ 1 - 7 | also observed in 



such operator A has (i) a finite index in Majorana repre- 
sentation (see Fig. [T]) - i.e. a finite number of occupied 
a-fermi states - or (ii) A has an infinite index (in the 
limit n — > 00). It was shown analytically [12j | for the 
transverse Ising chain (7=1) that in the case (i) the 
OSEE saturates in time while in the case (ii) the nu- 
merical results give firm evidence that the OSEE grows 
logarithmically as 



S = chit 



(4) 



where the coefficient c is the same for any infinite-index 
operator A = FB where F = w\ ■ ■ ■ w n and B is a finite- 
index operator. Thus we will eventually consider only 



the open noncquilibrium quantum XY chain 15( regard- 
less of the way the chain was coupled to the reservoirs. 
Therefore, the origin of such phenomen is attributed to 
the adjoint Hamiltonian H which gives exactly the uni- 
tary part in the master equations governing the de nsity 
matrix evolution of the boundary opened spin chain [14| . 

In the thermodynamic limit n — > 00 we expect that the 
dynamics will not change significantly if periodic bound- 
ary conditions in Majorana space are imposed. For finite- 
index initial operators being a product of Majorana op- 
erators near the center of the chain (Fig. Q]case (i)) the 
OSEE indeed agrees with the open boundary result which 
is also consistent with the area law as the area of the 
boundary remains the same. 
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FIG. 3: OSEE of |f for the XY(7,/i) model (HI with 
open boundaries (lines) and periodic boundary conditions in 
Majorana space (OSEE divided by two, symbols). Parameters 
are designated by The upper-most line denotes OSEE 

for the model with an additional term (|13[l with /1 = 0.4. 

For infinite-index operators however the area of the 
boundary doubles when periodic boundary conditions in 
Majorana space are imposed which results in a OSEE 
multiplied by a factor of two. Dividing the OSEE for 
periodic case accordingly, we obtain a perfect agreement 
with the results for the open boundary case as seen in 
Fig. [31 The upper-most line will be described later. 

The advantage of periodicity in Majorana space is that 
H can be diagonalized using the Fourier expansion 

followed by a Bogoliubov transformation. First we 
rewrite H in terms of u k which satisfy the same anti- 
commutation relations as Wj, 

n 

h = K"L"2fc-i + a fc4fc-i^2fe) ( 7 ) 
fe=i 

with a k — 2i(cos0fc — h) — 2"f sin cf>k- Introducing a± k = 
^(i-r^j-ua/s-i + u 2k) we diagonalize ([7|) and obtain a 
quasi-particle hamiltonian 

n 

H = J2 e k( & l^k ~ & -k & -k) ( 8 ) 

k=l 

with eigenvalues e k — |ofc| given by a dispersion relation 



e k =e(cj) k ) = 2^J{h-cos(j) k ) 2 + 1 2 sm 2 (j) k . (9) 

The OSEE is calculated by transforming u k back to Wj 
which gives us the time evolution of the canonical maps 
Wj{t) and therefore the evolution matrix <&, 

/ w 2 j-i(t)\ ( fl-j{t) 9l-j(t)\ fw 2 i-i\ nm 
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FIG. 4: Dispersion relation e(4>) (eq. [9} for XY model with 
7 = 0.5 and various h (indicated in the figure) and for the 
modified model (fl3|) (long-dashed curve). 

where fj , gj are real functions given as Fourier series 

Comparing the dispersion relation (Fig|4]) for (7, h) 
with the OSEE (Fig. [3]) we make the following 
Conjecture: The prefactor c in the logarithmic growth 
of OSEE for infinite-index operators and open bound- 
aries is given by the number m of stationary points of the 
quasiparticle dispersion relation e(4>), as c = m/6, if all 
stationary points are non- degenerate ]! aj '. 
We observe that the nontrivial stationary points where 
e'((j)) = only exist when \h\ < |1 — j 2 \ which determines 
the magnitude of the critical field separating the regions 
of m — 2 and m = 4, 

/ lc = |l- 7 2 |. (12) 

We verify the conjecture by considering more general 
Wigner-Jordan solvable spin chains, say by adding a term 
crJiT| +1 (7j_|_2 - or iu>2j-i u>2j+4 expressed with Majorana 
operators - to the XY- hamiltonian H in |T]), 

//' // • //»:<i<2- (13) 

3 

Choosing /1 = 0.4 we obtain a dispersion relation with 6 
stationary points (FigQ| which agrees with the result for 
the OSEE growth with the prefactor 6/6 (FigEJ. 

The idea of stationary points is supported by study- 
ing disordered chains where the Fourier transformation 
- which would result in the dispersion relation - is not 
applicable. We introduce disorder of strength e to either 
7 or h, as jj = 7 + Ej or hj = h + Sj, where Sj € [— e, e] 
are uniformly distributed random numbers. Although we 
shall only present data for /i-disorder, similar behavior 
has been observed also in the other case. 

We first check for a finite- index initial operator <t^ , 2 for 
which we find that the OSEE saturates in time, both in 
homogeneous and disordered cases (inset of Fig. [5]) . This 
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result is expected since for any finite-index initial prod- 
uct operator A (product of w/s) OSEE is upper-bounded 

by S(t) <K\og2 where K = Y%L i M\<"l™k\ A )$ is the 
Majorana index. The upper bound applies to any Hamil- 
tonian quadratic in Majorana operators and can be de- 
rived by a straight-forward generalization of the proof for 



the critical quantum transverse Ising chain in 12j . How- 
ever, no statement can be made on whether the OSEE for 
finite-index operators in the disordered model is higher 
or lower than in the corresponding non-disordered case. 

The OSEE of infinite-index initial operators is not 
bounded in general and the disorder significantly af- 
fects the dynamics of entanglement. Again, we will re- 
strict our interest to the simplest infinite-index operator 
F = w\W2 ■ ■ ■ w n as the results are qualitatively similar 
also for e.g. c^/ 2 or °"n/2' For a weak disorder, three 
stages can be identified in the evolution of infinite-index 
operators (Fig. [5]). Up to the time proportional to 1/e, 
the OSEE roughly follows the non-disordered case, after 
which it grows with a rate proportional to e until it fi- 
nally saturates to a plateau. The saturation phenomenon 
is a consequence of Anderson localization of eigenvectors 
of H in the Majorana space and must be contrasted with 
the logarithmic growth of the OSEE in time for the cor- 
responding non-disordered model. The plateau value of 
OSEE decreases with e although the quantitative relation 
cannot be established at present. 




FIG. 5: OSEE for the XY model with 7 = 0.2 in a randomly 
disordered magnetic field hj G [— e, e) averaged over 1000 real- 
izations for disorder strengths e £ {0.05, 0.2, 0.5, 0.8, 1.0, 1.5}. 
While main figure corresponds to the infinite-index operator 
F(t), the inset shows data for finite index operator 0^/ 2 (t). 



We have made an additional test where the external 
field is not randomly disordered but periodic such as 
hj = h + (— iyh' . Due to periodicity of hj we can 
again establish a dispersion relation which agrees with 
the logarithmic growth of the OSEE with the prefactor 
determined from the number of stationary points. 

In conclusion, we have demonstrated a discontinuous 



transition in the quantum XY-model at the critical field 
strength h c = 1 1 — 7 2 1 which separates two phases of dif- 
ferent entanglement production rate in the temporal dy- 
namics of operators. The transition exactly corresponds 
to a far-from-equilibrium phase transition in the steady 
state of the open XY chain [l~5[ . Beyond the XY model, a 
general relation has been conjectured which connects the 
prefactor of the logarithmic entanglement growth to the 
number of stationary points of the quasiparticle disper- 
sion relation. Interestingly, this temporal scaling of the 
operator space entanglement entropies is strongly remi- 
niscent of the size scaling of the ground state entangle- 
ment entropy found for the critical (gapless) XY models 
3 12 1 j where the prefactor c = m/6 is related to the 
number of sign changes m of the so-called symbol of an 
appropriate Toeplitz determinant. However, in our con- 
text the logarithmic growth (versus the saturation) of 
the entanglement entropy is not determined by gapless 
(gapped) nature of the Hamiltonian, but rather by the 
infinite (finite) index of the initial operators. 
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